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Abstract 

A detailed study of the relativistic classical and quantum mechanics of the massless harmonic 
^ ■ oscillator is presented. 
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I. INTRODUCTION 



One of the most important physical systems for both classical and quantum mechanics 
is the harmonic oscillator. In contrast to the nonrelativistic harmonic oscillator that is dis- 
cussed in most textbooks, the theory of relativistic harmonic oscillator is far from complete. 
The reason for this is the complexity of the problem related to the nonlinearity of differen- 
tial equations of motion for the classical relativistic oscillator. No wonder that even in the 
simple case of the massive one- dimensional relativistic oscillator there are problems with 

n 

identification of periodic solutions to equations of motion [1]. The problem of a quantum 
relativistic harmonic oscillator is usually formulated in one of three different frameworks: 
the Klein-Gordon, Dirac or Salpeter equations. The first one uses the spinless Klein-Gordon 
equation with a Lorentz invariant oscillatory potential j^J. However, the solutions of that 
equation are blamed by pathologies such as the appearance of ghost states. The second 
approach, referred to by Moshinsky [3] as the "Dirac oscillator" and describing spin one-half 
particles utilizes the Dirac equation with an appropriate combination of the scalar, vector 
and tensor couplings with an external field [4] . It can be successively applied to analysis of 
relativistic symmetries which recently were recognized experimentally in both nuclear and 
hadron spectroscopy |5j|. Unfortunately, this approach has no classical relativistic counter- 
part. Finally, the third approach follows from the relativistic Hamiltonian dynamics for a 
scalar particle and on the quantum level it is based on the spinless Salpeter equation jg]. 



The Salpeter equation is a "square root" of the Klein-Gordon equation [17] and can 

be regarded as its alternative . The serious advantages of the Salpeter scheme are the 
lack of problems with probabilistic interpretation on the quantum level as well as the classi- 
cally well-defined physical content of this theory. This last framework is frequently used as 
a phenomenological description of the quark- ant iquark-gluon system as a hadron model. 

Surprisingly, to our best knowledge, the simplest case of the massless relativistic harmonic 
oscillator was not discussed in the literature. In this work we perform a detailed analysis of 
the massless relativistic harmonic oscillator. In particular we find the exact solutions to the 
classical Hamilton equations as well as to the corresponding quantum Salpeter equation and 
discuss their basic properties. The article is organized as follows. In Sec. II, by integrating 
the corresponding Hamilton system we identify all kinds of possible motion of the oscillator 
as well as find its quantative characteristics. For an easy illustration of the dynamics of the 
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relativistic massless harmonic oscillator we also provide a graphical presentation of numerical 
integration of equations of motion. Section III is devoted to the quantization of the massless 
relativistic harmonic oscillator. 



II. THE ANALYSIS OF THE CLASSICAL RELATIVISTIC MASSLESS HAR- 
MONIC OSCILLATOR 

The Hamiltonian of the relativistic massless particle subject to the potential ^k 2 x 2 is 
given by 

H = c\p\ + ^K 2 x 2 , (2.1) 

where x and p are the position and the momentum of a particle, \p\ = \pp 1 ~ is the norm of 
the vector p (so c\p\ is the kinetic energy of the particle), k is a constant and c is the speed 
of light. Therefore, the Hamilton's equations are 

P 

X = C-. r, 

p = —K 2 X. 

(2.2) 

We point out that an immediate consequence of Eq. (2.2) is x 2 = c 2 , that is the length of 
velocity is c as should be for a massless particle. The familiar integrals of the motion in a 
central field [18| are the energy E and the angular momentum J: 

E = c\p\ + y* 2 , (2.3) 
J = x x p. (2.4) 

As a result of the conservation of the angular momentum J the motion is planar and we 
can restrict, without loss of generality, to the case of a particle moving in the (x 1 , x 2 ) plane. 
On passing to the polar coordinates x = (x l ,x 2 ) = (r cosy?, r sin <p) and p = (p l ,p 2 ) = 
(pcos8,psin6), where r = \x\, and p = \p\, we obtain from Eq. (2.2) the following system: 

r = ccos(9 — ip), 

tp = - sin(6* — ip), 
r 

p = —K 2 rcos(6 — (p), 

6 = K 2 -sm(9-(p). 
P 

(2.5) 



The integrals of the motion take the form 

E = cp + i/tV, (2.6) 
J = J 3 = rpsm(9 - cp). (2.7) 

From (2.5), (2.6) and (2.7) we find 
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r\jr 2 yE - —r 2 j - (Jc) 2 dip = ±\J\cdr. (2.8) 

We point out that the "+" and "— " signs correspond to the two possible orientations of 
the angular momentum. We choose, without loss of generality, the sign "+" and J > 
throughout this work. Now, from Eq. (2.8) we find that the trajectories should satisfy 

r (e - yr 2 ) > Jc, (2.9) 

and we can classify the types of motion as folows. The first two possibilities refer to J ^ 0. 
Namely, 



1) For r i^E — ^r 2 j = Jc, we get 
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r = J R=W— , (p = ut + (p , p = Po = — i = u t + (p o + -, (2.10) 

where u = % = ^-^ = c\/^k. So in this case we have a uniform motion in a circle with the 

R jn) y 2E 

linear speed \v\ = uR = c. This solution can also be obtained from Eq. (2.2) by demanding 
x-p = 0. Indeed, it can be easily checked that (2.2) and (2.3) imply the system 

^ j-i 3 9 2 

—x-p = E K X , 

dt 2 
— p 2 = -2k 2 x-p, 

d 2 2c 



— x = - — -x-p. 

dt \p\ 



(2.11) 



From (2.11) it follows easily that the orthogonality of x and p refers to the motion of a 
particle in a circle with radius \x\ — R — -y/ffl- 

2) For r i^E — ^r 2 ^ > Jc, the path lies entirely within the annulus bounded by the circles 
r = r min and r = r max , that is we have 

r m m < r < r max , (2.12) 



where r min and r max are the real positive solutions of the equation 

r 3 -tE + Jc = 0. (2.13) 



k2 j> 



2 

We find after some calculation 
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- - 2 V^ sin 3' ( 2 - 14 ) 



^(v^cosl - sin|), (2.15) 
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where sin a = ( |^ J , and < a < f . 



We now return to Eq. (2.8). An immediate consequence of integration of Eq. (2.8) is the 
the relation 



¥> = Vo + £ T , d% 2 =■ (2-16) 

- ''-I XJX(E-^X) 2 -(JC) 2 



In the particular case of r = r min 7^ 0, and r(t) > r min , t > 0, the integral from the 
right-hand side of (2.16) can be expressed by means of the elliptic integral of the third kind 
11(0, n, k) (see Ref. [19] . 3.137, Eq. 3), namely we have 



min V — min \ V max mm mm V mln 



(2.17) 

where r_ is the negative root of the polynomial from the left-hand side of Eq. (2.13) satisfying 

Train + ?" max + r_ = 0, T 2 _ > 7^ ax > 7^. (2.18) 

Clearly, Eq. (2.17) defines r as an implicit function of ip. 

Furthermore, for r = r max 7^ 0, and r(t) < r max , t > 0, the implicit equation for the 
trajectory can be obtained from (2.16) with the help of the elliptic functions of the third 

19[] , 3.137, Eq. 4). It follows that 



kind 11(0, n, k) and first kind F((/>, k) (see 

2Jc 



x II arcsin 



(r 2 - r 2 ) 

V — max/ 
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\ — min/V max / max min — / max min 



(iii 2 2 \ / 2 iv«2 ~\ jv>2 i "~ iy»2 ' \ / 2 
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h rl^F ( arcsin 4 / ^ o \ m ? , x I "T 2 m i " 

min 



\ V ( r max - r min)( r2 ~ ^) ' V T ~ ~ T 



To 



(2.19) 



Now, the four- dimensional system (2.2), where x = (x 1 ,^ 2 ), and p = (p l ,p 2 ) is completely 
integrable. Indeed, it possesses two integrals in involution E and J. Therefore, the motion 
between two circles with the radius r min and r max can be only quasiperiodic and periodic. 
Of course the case of the periodic motion refers to a closed path. This means that an angle 
A(p given by (see formula (2.17)) 

2 

Jc /" rmax dx 



n x\/x(E - \x) 2 - (Jc 
2Jc 




(2.20) 



should be a rational function of n, i.e. Ay? = nm/n, where m and n are integers. An 
example of a periodic motion between two circles is presented in Figs 1 and 2. It should be 
noted that the length of momentum (kinetic energy of the particle c\p\) has maximum at 
r = r min, decreases (increases) as r approaches r max (r min ), and for r = r max has minimum. 
Clearly, such behavior of the momentum of a massless particle is consistent with the form 
of (2.6). The values of r min and r max as well as extrema of the length of momentum can be 
expressed as a function of the energy by means of the implicit formulas (2.14) and (2.15). 
We finally remark that the case of the uniform motion in a circle discussed earlier [type 1) 
of the motion] refers to the condition r min = r max = R = \J~§^- 
The third type of the motion corresponds to J = 0, so we have 

3) r[E — -y-r 2 ) > 0. From this inequality we find < r < r max = ^H. On the other hand, 
taking into account (2.7) we find that for J = the system (2.5) reduces to 

r = ±c, 
= 0, 
p = =FK 2 r, 
= 0, 

(2.21) 

where 6 — ip = or \9 — ip\ = tt. Therefore a particle motion is uniform in a segment [0, -^^], 
more precisely, we have r = ±ct + r , where < r < and the two signs correspond to 
two possible directions of motion. Assuming that a particle moves in the x-coordinate line, 
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FIG. 1: The periodic solution of the system (2.2) obtained by numerical integration. The initial 
data are x = (0.479000,0.000000) m, p = (0.000000,1.290805) Jsm -1 , the parameter k 2 = 1 
Jm~ 2 , and c = 1 ms" 1 . 



(i.e. x = x 1 ), we get 

x = ±ct + x , -r max < x < r max , (2.22) 

where the turning points are x = r max and x = — r max . On setting xq = — r max we can write 
the trajectory explicitly as 

~2t 



x(t) = (-l)[*] c Jf_ (2 



1 1 !>■ 



(2.23) 



where T = is the period of oscillations of a massless particle between the turning 

points x = r max and x = — r max , and [a] is the biggest integer in a. The trajectory (2.23) is 




FIG. 2: The plot of the radius r = \x\ (solid line) and the length of momentum p = \p\ (dotted 
line) vs time obtained by numerical integration of (2.11). The initial condition is the same as in 
Fig. 1. 

illustrated in Fig. 3. Notice that at the turning points the momentum of a massless particle 
vanishes [see Eq. (2.6) for r = r max ] that is p min = 0. The maximum value of momentum 
Pmax — — is reached for x — 0. The time development of the momentum for p Q = and 
Xq = — r max can be written in the form 

p = (x 2 (t) - c^j , (2.24) 

where x(t) is given by (2.23). The plot of p versus t is shown in Fig. 3. Because the 
momentum of a massless particle tends to zero as its position approaches the turning point 
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FIG. 3: The plot of the coordinate (solid line) and the momentum (dotted line) of a massless 
oscillating particle vs time given by (2.23) and (2.24), respectively, where c = 1 ms _1 , k 2 = 1 
Jm" 2 , r max = 1 m, and T = 4 s. 

we deal with a "red shift" similar to the gravitational one. It should also be noted that the 
motion in the segment can be easily obtained from (2.11) by setting j||£L_ = ±1, that is x 
and p are parallel or antiparallel and therefore satisfy x x p = 0. Evidently, in the case of 
the system (2.2) this condition is equivalent to J = 0. We point out that the motion in a 
segment corresponds to the condition r min = and r max = for the nonnegative solutions 
to (2.13). We finally remark that the type of motion is completely determined by the values 
of the energy E and the angular momentum J. Namely, using the parametrization of r min 
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and r max denned by (2.14) and (2.15) we find 

3 3 

where j§ (llf) 2 = ^ re f ers to the motion in a circle, < j§ ( If) ^ < ^ corresponds to the 

3 

motion between two circles, and j§ (if) = ^' *' e ' ^ = ^ * s ^ ne condition for the motion in 
the segment. 

III. QUANTUM MECHANICS OF THE RELATIVISTIC MASSLESS HAR- 
MONIC OSCILLATOR 

In relativistic quantum mechanics the massless harmonic oscillator defined by the Hamil- 
tonian (2.1) is described by a massless version of the spinless Salpeter equation 

\h—i>{x,t) = ich^/^A x + —x 2 \ ip(x,t), (3.1) 

where A x = (^) 2 . Therefore the eigenvalue equation for the Hamiltonian H^e = E^ E 
takes the form of the pseudodifferential equation 



(3.2) 



Performing the Fourier transformation 



4,{x) = [ d 3 ke^^(k), (3.3) 

(27m) 2 J 

we get from (3.2) the following equation: 

2r \ 2F 
- A * + ^)' i ^ = ^MK), (3.4) 

where A k = (^) 2 - Finally, switching over to the spherical coordinates k = 
(k sin a cos f3, k sin a sin j3, k cos a), where k = \k\, and making the ansatz 

Mk) = X^-Yr(aJ), (3.5) 

where Y™{ot,0) are the spherical functions, we obtain the "radial equation" 

d 2 1(1 + 1) 2c ,\ 2E . . 



dk 2 k 2 (kH) 2 ) AV ' (kH) 
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To our best knowledge in the case of / 7^ the solution of (3.6) is not known. For I = the 
solution to (3.6) can be expressed by means of the Airy function Ai(x) 20], namely 



X {k) = CAi 



2c 



c 



(3.7) 



(2cKh)i 

where C is constant. We point out that / = was also the case discussed in [6j, where 
the recurrence was identified satisfied by coefficients of the formal power series expansion 
for the solution to the spinless Salpeter equation corresponding to the massive relativistic 
harmonic oscillator. Clearly, / = refers to the vanishing angular momentum, therefore we 
deal in this case with the quantization of the motion of a massless particle in the segment 
< r < r max = discussed in the previous section corresponding to the condition J = 
(third type of the motion). Furthermore, for I = the ansatz (3.5) takes the form 



47T k 



Demanding that ipE{k) is well defined for k = we find x(0) = (compare 21| Eq. (32.11)), 
which leads to Ai I I = 0. This quantization condition means that the values of the 



(3.9) 



(2creft)3 

energy E n , n — 1, 2, . . ., are given by zeros of the Airy function a n . We have 



E„ 



(2cnh) 



n 



1,2, 



Using the fact that the functions Ai(x + a n )/Ai'(a n ), n — 1, 2, . . ., where Ai'(x) designates 
the derivative of the Airy function Ai(x), form an orthonormal basis on the interval [0, 00) 
[l^ . we find that the normalized solutions (3.5) to (3.4) in the Hilbert space L 2 (R 3 ,ci 3 fc), 
with / = can be written as 




1 



-Ai 



2c 



-k + a n 



2vr(2c^)§ Ai(a n )fc \(2ckHY- 
From (3.10) and (3.3) we finally obtain the normalized wave functions such that 
V I I I 1 /' x „ . /,-/• . . / 2c 



(3.10) 



ll>n{x) 



11 ■ kv . . 
, . , dfcsm — Ai 

n vr (2ck?i)3 Ai (a n ) r J h \ (2ckH)\ 



■k + a. 



(3.11) 



where r = \x\. 

As in the case of the nonrelativistic harmonic oscillator with the probability density 
different from zero outside the turning points, the probability density p n (r) = \ip n (x)\ 2 
does not vanish for r > r mSLX (E n ), where r max (E n ) = ^ 2 ^ n , n = 1,2,...,. However, it 
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follows from the numerical calculation that p n (r) has no maxima for r > r max (E n ) (see Fig. 
4). Furthermore, taking into account all directions of the motion in the segment [0,r max ] 
classical limit does not deal with a single classical orbit but an ansamble of classical orbits 



22l | ) and taking into account that the probability of finding a particle in the spherical layer 



r, r + dr is inverse proportional to the surface of the sphere with radius r, we find that the 
classical probability density is given by the formula 



p c i(x) = p c \[r) 



9(r n 



47rr max r 



(3.12) 



where r max = and 6{x) is the Heaviside step function. Clearly, the normalization 
condition is of the form 



(fxp c i{x) 



p c \(r)dp(r) 



(3.13) 



where dp(r) = 4irr 2 dr. The comparison of the quantum probability density p n {r), and the 
classical one p c \(r) for r max (E n ) is shown in Fig. 4. As expected the differences between the 
quantum and the classical descriptions decrease as the quantum number n increases. 

We now discuss the expectation values of both the kinetic and potential energies. Using 



the identity 



23] 



1 f°° 
— — — - / xAi 2 (x + a n )dx 

i K)J 2 Jo 



(3.10) and (3.9) we get 



[A 



(i/j n \cpip n ) = c d 3 k\k\\ij n (k)\ 



-~a r , 



2 

~ —Er, 



(3.14) 



(3.15) 



where p = Hence, taking into account the form of the Hamiltonian in Eq. (3.1) we 

find 



A. 



(3.16) 



where f = We conclude that the virial theorem takes the nonstandard form in the 

case of the massless relativistic harmonic oscillator. More precisely, the roles of the kinetic 
energy and potential energies are exchanged. Interestingly, we have the same formulas on 
average kinetic and potential energies in the classical case. Indeed, from Eq. (3.12) it follows 
easily that 



cl 



y J r 2 p cl (r)dp(r) = ~E. 



(3.17) 



12 




).0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 
r [m] 




0.014 





FIG. 4: The plot of quantum probability density p n (r) = \if) n (x)\ 2 (solid line), where ip n (x) is 
the wave function (3.12) and the classical probability density p c \(r) (dashed line) given by (3.13), 



where r 



max — ' max 



{E n ) 



/2E„ 



. We set c = 1 ms 



-1 ,_2 



1 Jm 2 , and H = 1 Js. 



Therefore, by virtue of the first equation of Eq. (2.3) we have 

2 

(cp)ci = ^E. 

We finally write down the following approximate relation obtained numerically: 

' max 

(E n ) \/2E n 



2k ' 



(3.18) 



(3.19) 
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where the formula is exact in the limit n — > oo. The approximation in (3.19) is very good. 
The maximal relative error \((i/) n \ril) n ) — r max (£' n )/2) / (i/j n \ri/j n ) | arising in the case with n = 1 
is about 5% and is lesser than 1% for n = 2. The fact that the limit n — > oo when we have 
the exact equality in (3.19), is the classical limit is confirmed by the classical formula 

rp cl (r)dfi(r) = = ^ (3-20) 

following directly from Eq. (3.12). 



IV. CONCLUSION 



In this work we study the relativistic massless harmonic oscillator in both classical and 
quantum cases. It seems that the obtained results concerning the classical oscillator are of 
importance not only from the physical point of view. Indeed, Eq. (2.2) is one of the simplest 
examples of a nonlinear Hamiltonian system with constant length of velocity. As far as 
we are aware such an interesting class of nonlinear dynamical systems was not discussed in 
the literature. Referring to the observations of this work related to the quantum mechanics 
of the relativistic massless harmonic oscillator we wish to point out that Eq. (3.11) is, 
to our best knowledge, the first example of the nontrivial exact solution to the Salpeter 
equation. We also stress the good behavior of the corresponding probability density and 
expectation values of observables which confirms the correctness of the quantization based 
on the massless spinless Salpeter equation. Furthermore, we obtain the exact formula (3.9) 
on the spectrum of the Hamiltonian. It should be noted that for the Salpeter equation only 
energy bounds were analyzed in the literature so far (for the massive relativistic harmonic 



oscillator see Ref. 



10|). Finally, we have obtained the interesting form of the virial theorem 



for the massless relativistic harmonic oscillator with the exchanged roles of the kinetic and 
potential energies. 
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